Commuting Matrix Solutions of PQCD Evolution Equations 



Mchrdad Goshtasbpour 1 and Seyed Ali Shafici 2 

1 Dept. of Physics, Shahid Beheshti University, G. C, Evin 19834, Tehran, Iran. 
2 Biophysics- Biochemistry Department, Rafsanjan University of Medical Sciences (RUMS), Rafsanjan, Iran. 



CO 

o 

CN 



43 

i , 

Or 
<D ■ 
43 



> 

in 

oo 

co 
co 

o 



x 



A method of obtaining parton distributions directly from data is revealed in this series. In the 
process, the first step would be developing appropriate matrix solutions of the evolution equations 
in x space. A division into commuting and non-commuting matrix solutions has been made. Here, 
well-developed commuting matrix solutions are presented. Results for finite LO evolution match 
those of standard LO sets. There is a real potential of doing non-parametric data analysis. 



PACS numbers: 12.38.Bx 



INTRODUCTION 



Given conditions on parton distribution functions 
(pdfs), coming from existing measurements in the accel- 
erators, e.g. from a measurement of the structure func- 
tion F 2 in DIS, through a fixed order (LO, NLO, ...) 
formula from PQCD for F 2 , or a given cross section for 
pp scattering in a similar manner, can be considered the 
initial conditions for the pdfs in the intgro-differcntial 
DGLAP evolution equations of the same order in per- 
turbation; based on which solving the DGLAP equations 
generally allows PQCD prediction for processes of inter- 
est in accelerators, and eventually in the colliders. Thus, 
the method of solution of DGLAP proposed in the year 
2000 and beginning to be presented in details here,* 
can be of serious interest. For simplicity the method is 
illustrated at LO PQCD. 

DGLAP evolution equation is of the following simple 
form for the non-singlet quark distribution Jns'- 



df NS (x,t) 
dt 



f ^P qq (-)f NS (y,t) 
Jx v y 



(i) 



Usually, the derivative in ([1} is with respect to Ln(Q 2 ), 
where 



dt 



dLn(Q 2 ) 



a s (Q 2 
2ir 



(2) 



In ©, and throughout this work, the LO value of the 
strong coupling constant a s (Q 2 ) is used, either in the 
historical approach based on Aq C£) , or in the present day 
approach based on experimental value of a s {M z ) directly, 
[a, l6( , equivalently (with no difference) . 



n DIMENSIONAL SPACE, PDF VECTORS, 
SPLITTING FUNCTIONS MATRICES 



for large x% < 1. A set of n basis vectors for each pdf 
set can be defined on each x%. The splitting function, 
in terms of ([3]), can be calculated as an n x n matrix. 
The choice of x\ and n in ([3]) is eventually constrained 
by a close fit to the x points of the data, for better 
interpolation. 



Convolution integrals 



Convolution integrals are of the form >t 



P<8)/: 



in which the simplest parton distribution function 
f(x,t) = fNs(x,t) in the discrete basis is an n-tuple, 
whose ith element is f{xi,t). As the general convolution 
integral is made a discrete sum, P^,i,j = q,g within 
the kernel P, see (fT4|). become lower triangular (l.t.), 
due to the limits of integration of the resulting sum, and 
banded (bdd), due to divisions of Xi in sec (0 to 
(|10p . The resulting algebraic simplicity, i.e, commutativ- 
ity, of l.t. bdd matrices is the essential element of the 
following analysis of the solutions of DGLAP equations. 



A hybrid finite difference computation of splitting 
functions, an example 

To calculate the matrix form of a splitting function, 
e.g. P qq in (pj, a combination of discrete finite difference 
approximation of the unknown pdfs and continuous 
integration of the known splitting functions is mixed 
in our method of evaluation of the convolution, as 
follows. Then, (P qq )ik is extracted from the coefficients 



of fk = f(x k ,t). 
equivalent forms 



At LO, P qq is given in the following 



Define discrete Bjorkcn x%, 

Xi = (xi)V = I: -, n , 



(3) 



Pa. 
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[i + y 2 " 
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f l + y 2 
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[ 1-2/ . 


+ ~3 


[(1-2/)+ 



(5) 
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Given the definition of " +" rcgularization, aimed at re- 
moval of the infinity at y = 1 in the kernel P qq in the 
context of the convolution integral, 



dxf(x) 



1 - x 



(6) 



convolution in (|T|) becomes: 



i !/ y3Vi-y/+ Jx v 



4/1 
3 



(7) 



= h+h- 



In the left integrand of (O, an extra y is placed to have 
the more useful parton momentum distribution function. 
Note that the DGLAP equation can be written for mo- 
mentum distribution; i.e., in it f(x) — > xf(x), by putting 

an extra x on both sides, and a factor of v.- = 1 within 

' a y 

the integral. Different techniques of integration may be 
used in DGLAP. Here, we use only integration by parts: 



« = /(-) 

y 



#(f) 

dw = —dy, 

dy 



4 1 + y 2 /" 

cfo = -(- )dy ^>v= / dv; 

3 1-2/ 7 



7x = /(x)w(l) - /(l)w(a:) 



i <2/(|) 
v{y)—f ! -dy 



dy 



(8) 



1 ,x <ff(y) 
v{-)—r-dy, 



f(x)v(l) , 

/ 2 = -/(x)[«(l)-«(0)] 

/ = /(xH0)+ [\(*)VM_dy. 

Jx y dy 



Here and later, when convenient, the variable change 
^ — > y is used. Bringing in a finite difference approxi- 
mation of the differentials: 



i k ~ 1 

I = v(0)f( Xl ) + J2 f v{xi/y) 

k=l ~, 



fk — fk-1 
Xk - Xk-1 



dy 



Xi-1 



>{xi/y) 



fi - ft- 



Xi X%—\ 
fk — fk-1 



— dy 



Y] [ v{xi/y)— — -^-dy 
f— f J %k- Xk-i 

k—1 „ 



v(0)f(xi) 



x k - x k ^i 
v(xi/y)dy 



(9) 



fi-i 



Xi Xi — i 



/ v(xi/y)dy + Y] [ v ( x i/y)— — ^^-dy. 
J J Xk ~ Xk ~ l 

Now, a bdd l.t. P qq can be read (extracted) from the 
coefficients of fk ■ 



(P qq )u = «(0) + 



1 



P qq -< 



{Pqq)ik 



Xi Xi — 1 
xi 

1 



X k - X k -1 



v(xi/y)dy, 



v( Xi /y)dy (10) 



Xk 



Xk+1 ~ X k 



v(xi/y)dy. 



Xk + l 



The other three kernels are essentially obtained in a 
similar manner and presented in Table HI whose third 
column contains a small example set of numerical results 
whose x— points set is of actual use for NMC [§]. 



COMMUTING MATRIX SOLUTION OF DGLAP 
IN X SPACE FOR FINITE Q 2 INTERVAL 

The present section contains details of analysis of com- 
muting matrix solutions of DGLAP equations ([Lj) and 
(TH)) , for finite Q 2 interval. 



Non-singlet 

Using discrete form of the splitting function (flU)) , dis- 
crete form of ([1]) is: 



dt 



£P y /f s (t)^/(t) = P./(t). (11) 
i=i 



3 



Splitting Function 


commuting matrices 


Numerical 


p - ifkkvl), 

r 1t — 3 V- 1-y > + 


(P qq )a ~ v qq (0) + x% _ Xi _ 1 J v qq (xi/y)dy 

x i 

x k-l x k 

(P qq )ik ~ Xk -l k _ 1 f v qq (xi/y)dy Xk+ \_ Xk J v qq (xi/y)dy 
x k x k+i 


{-3.96,2.54, 
0.62,0.29, 
0.16,0.10, 
0.06,0.04} 


A 9 = |[y 2 + (i-y 2 )] 


x i-i 

(P qg )u - v w (l) + Xi _ Xi l J v qg (xi/y)dy 

x k-l x k 

{P qg )ik - Xk -l k _ 1 I v gg (xi/y)dy Xk+ l_ Xk J v qg ( Xl /y)dy 

x k x k+l 


{0.06, 0.07, 
0.04,0.03, 
0.02,0.02, 
0.01,0.01} 


p _ 4rl + (l-i/) 2 1 
r m 3 I y -1 


(Pg q )u - v gq (l) + J v gq (xi/y)dy 

x i 

"h-1 Zfe 

(Pgq)ik - Xk J Xk _ 1 I v aq {x t /y)dy Xk+ l_ Xk f v gq {x l /y)dy 
x k x k+i 


{0.24, 0.49, 
0.56,0.63, 
0.70,0.75, 
0.79,0.82} 


^ + S(l-y)( 3 -^) 


x i—l 2 

{Pg 3 )n ~ a gg + Xt -l l _ 1 I v gg (x l /y)dy, a gg - g 1 

x i 

x k-l x k 

(Pgg)ik ~ Xk - Xk _ 1 I vgg(x 1 /y)dy [_ Xk f v gg (x l /y)dy 
x k x k+i 


{-9.43,6.38, 
2.41, 1.90, 
1.80, 1.81, 
1.84,1.88} 



TABLE I. column 1: LO splitting functions, column 2: Commuting splitting function matrices, with Vij = f Pijdy, for ij = qg 

or gq, v qq = v in (j8j), and v gg = + 6< ' 1 ~ y ^ + 6j/(l — y). Column 3: first column of the l.t., bdd splitting function matrices 
calculated from column 2, with a choice of n — 8, and xi = .712. 



Independence of the kernel P = P qq from the variable 
t leads to the following solution to the matrix equation 
(ITTj) for finite evolution from to to t 

f(t) = e^ p f(t Q ) ee P t f(t Q ). (12) 

To obtain the last equality in (fl"2j) . the diagonal and non- 
diagonal elements of the l.t. matrix P can be separated 
into two commuting matrices as P = Pa + Pol, with 
obvious definitions, including po = (P qq )u. Thus, being 
strictly l.t., Pa" — , and this in turn implies a finite 
expansion of the exponential as: 

p t (i-i ) = e (t - to)Po /.E^- ( 13 ) 

Therefore, the finite evolution machinery, for the non- 
singlet solution of DGLAP, is constructed with the final 
(t — to) dependent, ready for computation, P t matrix. 



Singlet 

DGLAP evolution equation has the form of two 
coupled equations for the singlet quark distribution S, 
coupled to gluons: 



d_(E(x,t)\ 
dt \g(x,t)J ~ 

(14) 

(|14p can be cast into a matrix equation as: 

d -^= f 1( *Ps£)f(y,t). (15) 
9t Jx v y 

Operating on the 2n dimensional external product 
space of the coupled pdfs / , discrete form of splitting 
functions are placed in a (2n) x (2n) matrix: 

P s = I Pqq 2UfPq9 ) . (16) 
\ P gq Pgg J 

Based on arguments similar to those for ([TTj) . (JTTJ) gives 
a construction of the matrix equation of (|15[) . leading to 
its finite (t — to) dependent solution, as in the case of 
(fl"2"|) . due to independence of the kernel from t. 

f s (t) = P s .f S (t) => f(t) = e ( *- to)Ps / S (io). (17) 

The kernel Ps is expressed in an external (Kronecker or 
direct) product, (g), space of a 2 x 2 matrix space and an 
n x n matrix space of bdd, l.t., thus commuting, splitting 
functions. Ps is expanded, as a direct product, in terms 



4 



of the four 2x2 identity and Pauli matrices: 




as follows: 



Ps = I ® Po + P, 



Po 



(Pqq + Pgg): 



i 1 (18) 

P = O3 (8 P 3 + -02 ® (2n/P gs - Pgq) + -<7i (8 

(2n/P gg + P 3g ) = ct 3 (8) P 3 + 2n/o- + ® P 9g + <j_ <8 P gq , 



P3 o qq ^ 9g ^' 



where a± = \{o~\ ± £02)- Thus, 



,{t-t Q )Ps 



,(t-t )I®Po Jt-t )P 



(19) 



(fT9|) holds as the commutator [Pol, P] = 0, as an in- 
stance of easily proved [7 £g> a, A <g> 6] =0 for the direct 
product of m x m identity matrix 7 and general matrix 
A with any two bdd lower (or upper) triangular, thus 
commuting, n x n matrices a and b. 

Equation (|19j) may be written in its matrix form, ready 
for computation: 



e (t-t )Ps = e 



71=0 



/ ® e ( *- io)p ° [7 ® Cosh((i - t )P)- 
7®S , m/t((t-i )P)P- 1 P]. 



(20) 



In ([20]), P = VP 2 , is defining the length of P, like a 
2-norm, except that its components (vector coordinates) 
are not numbers but commuting bdd, l.t. matrices made 
of splitting functions. P 2 = P 2 is a system of equations 
to be solved for computing P. The last term in (|2TJ|) can 
be put in an alternative form: Sinh((t — to)P). Separat- 
ing the diagonal and the strictly l.t. part of the splitting 
functions Po and P in ([2TJ]) , brings the possibility of hav- 
ing a finite sum for the analytical evaluation of the singlet 
during the computation, similar to what was witnessed 
for the non-singlet, (|T3)) . 

Equations (Q~3|) and ([20|) are called analytical commut- 
ing matrix solutions as they are analytically derived com- 
muting matrices (finite sum representations of exponen- 
tials of commuting bdd l.t. matrices - computed nu- 
merically) for finite Q 2 interval solutions. In contrast, 
to bypass interpolation in data analysis in jij, a non- 
commuting matrix solution is formulated, using exist- 
ing numerical routines for computation of exponentials 
of non-banded, thus non-commuting, triangular matri- 
ces. 



COMPARISON OF PDF RESULTS FOR FINITE 
EVOLUTION OF NS-S-g WITH MSTW AND 
DISCUSSION 

The machinery thus far has equipped us with the ca- 
pability to take a discrete set the initial or input pdfs, 
at x-points patterned according to ((3]), at to from some 
present standard LO set; evolve it ourselves to t (our out- 
put), and then compare our results with their (final) pdfs 
at t. This is actually done via a quick computation, FIG. 
[1] [2] with a small n = 180, xi = .95, (|3|), input taken at 
Ql = lOGev 2 from MSTW, fl. 

Our pdfs are assumed to be everywhere in the limit 
of high energy zero quark mass ST_/p(5) symmetry. Five 
flavor non-singlets and singlet, shown in FIG. [TJ [2J have 
definitions: 



qz = u — d, qs = u + d — 2s, 

Q15 = u + d + s — 3c, <724 = u + d + s + c — 46. 

q s = S = u + d + s + c + b; 

q = qtot + q, q = u,--- ,b. 



(21) 



In addition, three non-singlet valence quarks are as- 
sumed along with MSTW: 



qv = qtot 



<1 



s: 



(22) 



and are directly evolved; furthermore, all the sea quarks 
evolution is inferred from ([2~T1 12"21 ) for the figures. 

There are two further assumptions. It is essential 
to have a continuous coupling constant, a s , as flavor 
number changes, as Q 2 crosses the quark mass m 2 = 
22.56 Gev 2 , [6(, in the finite interval of evolution, SQ 2 = 
[10,1000] Gev 2 . The LO value of a s (m z ) = .13939 of 
MSTW is used here together with the related considera- 
tions Absence of this continuity constraint, brings 
inconsistency and considerably increased deviations. We 
also assume a simple continuity of pdfs as the flavor num- 
ber changes at the same Q 2 = m 2 . We are not claiming 
that there is no discrepancy between our assumptions 
and those of MSTW. 

The quality of match of our pdf outputs at Q 2 = 
lOOOGev 2 (magenta dots in comparison) with MSTW 
(solid black lines) can be observed in FIG. HJ [2] Similar 
matching with any other good solution of DGLAP equa- 
tions at LO is expected. We have had explicit trials with 
CTEQ5L (I and GRV 0] pdfs. InFIG.H three indirectly 
evolved pdfs: u, d, s, are not compared to MSTW. Their 
absolute deviations from MSTW pdfs are very small and 
very similar to that of the c quark, while very different 
from the 6 quark, FIG. [5l |H1 Thus, their comparison is 
just as good as that of the c quark. Close scrutiny of 
FIG. [2] shows the difference of fit between the two, c and 
6 quarks. There will be a related discussion of b problem 
in the ending paragraphs. 

Quantitative description of the fit , depending on choice 
of the beginning point x\, may be described, at each 



5 



point x, (j3j, in terms of absolute deviation, S, defined 
as deviation of our result from the corresponding (final) 
MSTW pdf; or relative deviation, Sr, denned as 5 divided 
by the value of the (final) MSTW pdf. It can be seen that 
our method based on discrete x, ©, induces oscillations 
in the deviations. 





1.4 : 




1.2 1 
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MSTW 1000 Gev 2 




our output 
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from top large x: _ 




<7i5, <lg, U v , Ov 


>A 0.4 '- 
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-4 -3 -2 -1 

FIG. 1. Comparison of six evolved valence and valence 

dominated non-singlet pdfs with MSTW, vs. log\o(x), at 
lOOOGeu 2 . 




-4 -3 -2 -1 



FIG. 2. Comparison of gluon and four sea or sea influenced 
evolved pdfs: b, c, singlet, and 924 with MSTW, vs. log\o(x), 
at lOOOGei) 2 . Three evolved pdfs: u,d,s, axe not compared 
due to triviality, see FIG. [S] and explanation in the text. 

The order of calculation of evolved pdfs is according 
to increasing i, or decreasing xi, in ([3]). As Xi decreases, 
three areas may be differentiated. Relative deviations for 
large x area are very large, which is not far from expec- 
tation due to very small size of denominator pdfs there, 
and having largest finite difference intervals in x, with- 
out smoothening effects of summation, of the convolution 
integral, which is just beginning there. 

The second middle to small x area approximately cor- 
responds to x € (1 0~ 4 , VlO) (stricly speaking, for near 
zero gluon and s v quark there VT0 — > FIG. [3]). 

In this area, for u v ,d v ,qa,qi5 which have the smallest 




FIG. 3. Relative deviations of all directly evolved nine NS- 
S-g pdfs with respect to MSTW. 




FIG. 4. Relative deviations of four directly evolved pdfs with 
smallest, and one, 524, with largest maximum Sr for x £ 
(10~ 4 , VlO) (see text for VW)). 

relative deviations, FIG. [3j maximum of 5p. FIG. [J is 
below .0007 (neglecting the last small x oscillation of d v ). 
For s v , FIG. [3j Sr is somewhat larger and problematic, 
due to smallness of s v , and that both input and output of 
MSTW cross zero in the middle of the second area. For 
924, 5r is the largest, FIG. [3J by a considerable margin; 
its maximum is below. 007, FIG.|U at x « 10 -4 . It may 
be expected from application of finite difference method, 
that deviations and maximum Sr decrease with increase 
in x\. We found this to be the case in a limited trial of 
letting x\ = .95 to take powers of 1/4, 1, 2. 

The exceptional behavior of (724 is related to different 
behavior of b quark from the other sea quarks, FIG. [3151 
[7] Further investigation is demanded; it may be due to 
difference of our assumptions and MSTWs, e.g., it may 
be related to flavor number scheme. 

In the third area of deeper small x, approximately x < 
10~ 4 , Sr ceases to be as good. We may return to details 
of the subject in future. 

In conclusion, we have developed a detailed numerical 
presentation of a commuting matrix solution of DGLAP 
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FIG. 5. Absolute deviations of the five indirectly evolved sea 
pdfs of FIG. [2 from MSTW. 
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FIG. 6. Absolute deviations of five evolved pdfs with respect 
to MSTW. It indicates a relation between the problems of 
large positive 6b and large negative Sq24, posed in the text. 



evolution equations at LO PQCD, as a first step of realiz- 
ing the potential of doing non-parametric data analysis. 
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